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ABSTRACT 



Context. Converging flows with strong magnetic fields of different polarity can accelerate particles through magnetic 
reconnection. If the particle mean free path is longer than the reconnection layer is thick, but much shorter than the 
entire reconnection structure, the particle will mostly interact with the incoming flows potentially with a very low 
escape probability. 

Aims. We explore, in general and also in some specific scenarios, the possibility of particles to be accelerated in a 
magnetic reconnection layer by interacting only with the incoming flows. 

Methods. We characterize converging flows that undergo magnetic reconnection, and derive analytical estimates for 
the particle energy distribution, acceleration rate, and maximum energies achievable in these flows. We also discuss a 
scenario, based on jets dominated by magnetic fields of changing polarity, in which this mechanism may operate. 
Results. The proposed acceleration mechanism operates if the reconnection layer is much thinner than its transversal 
characteristic size, and the magnetic field has a disordered component. Synchrotron losses may prevent electrons from 
entering in this acceleration regime. The acceleration rate should be faster, and the energy distribution of particles 
harder than in standard diffusive shock acceleration. The interaction of obstacles with the innermost region of jets in 
active galactic nuclei and microquasars may be suitable sites for particle acceleration in converging flows. 
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1. Introduction 

The interaction of outflows with themselves or their 
environment usua lly leads to particle acceleration (e.g. 
iRieger et al.l l2007f) . and the most common acceleration 
mechanism in these circumstances is tho ught to be dif- 
fusive sho ck (Fermi I) acceleration (e.g. Krvmskii Il977t 
[Be!il[T978allbl : iBlandford k Ostrikerll978l:lDrurvlll983ri . This 
process works well when a plasma with a weak disordered 
magnetic field (B) suffers a strong shock. The Fermi I accel- 
eration efficiency depends on the energy gain when particles 
cross the shock, and on the escape probability downstream 
the shock, which is related to the shock compression ratio 
R = v/v' , v'(v) being the post(pre)-shock velocity. 

In strong non-relativistic adiabatic shocks in the test- 
particle approximation R equals 4, and the res ulting energy 
distribution is Q(E) cx l/E 2 fe.g. lDrurv|[l983h . In radiative 
shocks, R 3> 4, but densities become very high and energy 
losses and damping of magnetic irregular ities may limit the 
efficiency of the acc eleration process (e.g. lDrurv et al.l[l996L 
iReville et al.ll2007l ; see nevertheless Sect. 12. 2p . For plasmas 
that carry a strong perpendicular B-component, R is ~ 1, 
and Fermi I acceleration is suppressed. However, if the B- 
field has reversals, energy can be dissipated via magnetic 
reconnection in regions where the charge density is not 
sufficiently high to sustain different polarity _B-lines (e.g. 



Parkerlll957t iPetscheld 119641 ISpeiserlll965t ISon nerud 119711. 
Zenitani k Hoshinof 2001 1. allowing for R ^> 1 (e.g. iDrurvl 



2012). 



Studies show that magnetic reconnection can accel- 
erate the bulk of particles up to the pre-reconnection 
Alfvenic sp eed, with some particles reaching even higher en- 
ergie s (e.g. IZenitani k Hqshinoll200ll : iLvubarskv k LivertsI 
2008; iKowal et al.l 1201 lh . Magnetic reconnection may 
also lead to the formation of strong shocks, accelerat- 
ing particles through the standard Fermi I mechanism 
in an otherwise weakh 
Il98l iBlackman k Field! 

Fermi I particle accelera tion may operate in the 



pressive flow (e.g. 1 


^orbesl 


iSironi k Spitkovskv 


20111). 



recon- 



necti o n region i t self ( e . g. Ide Gouveia dal Pino k Lazarian] 



120051: iGianniosI 120101: (DruryJ |2012j), since particles can 
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bounce back and forth between the layer and the converging 
flows. The mean free path A of these particles can eventu- 
ally overcome the layer thickness, AY, and then particles 
will not interact significantly with the reconnection layer, 
but with the pe rpendicular B-lines of the converging flows 
(|Gianniosll2010f) . In this case, the advection escape proba- 
bility becomes formally zero because the flows do not move 
away from the reconnection layer. This is true as long as 
the mean free path of the particle parallel to B, A|| , is much 
shorter than the characteristic size of the entire reconnec- 
tion region, AX (to order the magnetic field reversal scale) . 
Escape in the opposite direction to the incoming flow can 
be neglected unless the region extension in this direction is 
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Fig. 1. Sketch of the considered scenario: a) the interac- 
tion of one flow carrying different polarity i?-lines with an 
obstacle (see Sect. [3]); b) the interaction of two flows with 
different polarity B-lines. 



<C AX. A schematic picture of this scenario is presented in 
Fig.Q] which shows two cases: the interaction of a flow car- 
rying different polarity £?-lines with an obstacle (described 
in Sect. [3]); and the interaction of two flows with different 
polarity i?-lines. As shown in the figure, particles can spiral 
back and forth between the two sides of the reconnection 
layer until they diffuse away from the region. 

In this note, we analytically explore the Fermi I mech- 
anism when particles interact with converging flows in the 
null-escape probability regime. We assume that particles 
can effectively diffuse through the incoming flows. We fo- 
cus on the case when magnetic reconnection is the mech- 
anism that leads to R ^> I. We derive in Sect, [^analyt- 
ical estimates for the acceleration rate, the maximum en- 
ergy and energy distribution of the accelerated particles, 
and present in Sect. [3] an illustrative case in the context 
of galactic and extragalactic jets. For simplicity, the dis- 
cussions are restricted to Newtonian flows, although our 
conclusions qualitatively apply to the relativistic case as 
well. For some magnitudes, we will adopt the convention 
Qx = (Q/10 x cgs). 



2. Particle acceleration in converging flows 

In converging flows, particles slowly diffuse perpendicularly 
to the J5-lines, with Aj^ ~ r g /x, where x, with value > 1, 
is the total to the disordered magnetic energy density ra- 
tio. Remarkably, a disordered B-componcnt can enhance 
the reconnection rate an d the reconnecting flow veloc ity v, 
as shown for instance in lLazarian fc Vishniad (|1999f) . The 
following limit for Fermi I to operate between the converg- 
ing flows can be imposed: in the current sheet, A^ > AY, 
i.e. -Emm ~ \qB r AY, where B r is the B-field in there. 
Although one can expect that B r <C B close to the center 
of the current sheet, we conservatively assume that over the 
whole reconnection layer B r ~ B. On the other hand, the 
_B-lines become more and more entangled close to the re- 



connection layer, so one can set there x = !• The previous 
condition becomes 

E> £ min ~0.3B AF 9 TeV. (1) 

In addition, to be efficiently accelerated, particles in the 
incoming flows should not drift too early along the B-lines, 
parallel to the layer, out of the reconnection region (cross- 
field diffusive escape will be slower by 1/x 2 ). This is a sort 
of Hillas limit: A|| ~ x r g < AX, or 

E^ ax ~qBAX/x^3B AX loX - 1 TeV. (2) 

This gives the maximum dynamical range 

E^jE min ~ 10 x" 1 (AX 10 /AY 9 ) . (3) 

iGianniosI (|2010fl accounted for the electric field between 
the sides of the reconnection layer: f ~ «B. The e-field 
accelerates the particles perpendicularly to B and v, with 
£ max ~ (v/c)qBAX. The fraction of time spent by par- 
ticles within the reconnection layer during an acceleration 
cycle is AY/ct cyc \ c , and thus this effect will be important 
at E ~ i? m in or for very ordered fields. 

From Eqs. (J5H3J) one can see that the quantities X an d 
AX/ AY are crucial for the efficiency of the discussed accel- 
eration mechanism. The derivation of X requires a detailed 
treatment of the plasma, but as mentioned, given the dra- 
matic change of the -B-geometry, close to the reconnection 
layer % ~ 1. On scales of AX, and given the high plasma 
characteristic speeds, the disordered i?-field generated in 
the reconnection layer may propagate to the converging 
flows, although such a ^-component may also have an ex- 
ternal origin. To estimate AX/ AY, we briefly discuss the 
flow dynamics below. 

2.1. Flow dynamics 

The plasma that enters, flows along, and eventually leaves 
the reconnection layer, as shown in Fig. [51 can be described 
by the particle and energy flux conservation equations: 

Ap v = A'p' v' and (4) 

AvB 2 /8ir = A'p'v' 3 /2, (5) 

where unprimed/primed quantities correspond to the in- 
coming/outgoing flows, p is the density, v the velocity, and 
A ~ AX 2 and A' ~ AAXAY are the in- and out-flow ar- 
eas, respectively. The flow pressure is assumed to be negligi- 
ble in the converging flows because of strong B-dominance. 
When leaving the region, under negligible external pressure 
and reasonable geometries, the flow kinetic energy becomes 
dominant and pressure can be neglected as well. 

Given the symmetry of the problem, there is no momen- 
tum conservation equation, and the conditions deep within 
the reconnection layer are hard to determine. Equations 
S] and \5\ yield v' « v A — V 1 B 2 /inp and thus p' « 
p(A/A')(v/va), but one has to make additional assump- 
tions to obtain AX/ AY. An estimate may be derived as- 
suming one-dimensional adiabatic compression of B r , the 
remaining field after reconnection in the layer, until it bal- 
ances the incoming field pressure. This gives p' / p ~ (B/B r ) 
and therefore an area ratio A' J A = (B t /B)(v/va), or 
AX/ AY - 4{B/B T )(v A /v). An upper limit for AX/ AY 
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Fig. 2. Sketch of the converging and the outgoing flows in 
the reconnection region. 



would come from AY ~< r g >, where < r g > is the gyro- 
radius of the average particle in the reconnection layer. A 
lower limit may be derived assuming equilibrium between 
the upstream magnetic pressure and the thermal pressure 
inside the layer of thickness AY, implying (AX /AY) ~ 
A(va/v). One can therefore conclude that the dynamical 
range may span several orders of magnitude in energy, 
but it needs an efficient magnetic-to- kineti c energy tran sfer 
(which may be reasonable, as suggested in lDrurvl[20rl . 

2.2. Acceleration and radiation processes 

Like in standard diffusive shock acceleration, particles gain 
energy by moving back and forth between the incoming 
flows as AE / E ~ v/c per cycle. In addition, since there 
is no advection directed outwards, the probability to cross 
back to the other side of the reconnection layer becomes 1 
after covering few in each incoming flow, i.e. £ C ycic ~ 
(k/x) f g /c 7 with k ~ 10. This i C ycie gives an acceleration 
rate 



E n 



AE/t cyclc = ix/k) [v/c) q Be, 



(G) 



and the zero advection escape probability renders a parti- 
cle energy distribution Q(E) oc 1/E up to energies a few 
times lower than E^ x , at which the distribution drops very 
quickly. Note that this mechanism is faster in E acc ((v/c) vs 
(v/c) 2 ) and yields a harder Q(E) (1/E vs 1/E 2 ) than stan- 
dard Fermi I. Such a hard distribution of accelerated parti- 
cles can become dominant in pressure. This should smooth 
the v profile a nd affect Q(E ) (as in non-linear Fermi I ac- 
celeration; e.g. lDrurvlll983l ). 

Radiation losses can stop the acceleration at E < 
-^max- Typically, for protons these losses are fairly in- 
efficient, whereas for elect rons synchrotron cooling (e.g. 
Blu menthal fc Gouldlll970l) tends to be important, which 
renders a maximum energy 



E sy 



60 ( X /k) 1/2 (v/c) 1/2 B- 1/2 TeV. 



(7) 



In environments with very dense matter or photon fields, 
other cooling processes might be relevant, such as rel- 
ativistic Bremsstrahlung and inverse Compton (IC) for 



electrons (e.g. iBlumenthal fe; Gouldl Il970l ). and proton- 
proton collisions (pp), synchrotron and photomeson produc- 
tion for protons (e.g. IKelner et all 120061 : lAharonianl 120001 : 
iKelner fc Aharonianl 120081) . Diffusive escape parallel to B 
can also be more restrictive than the Hillas limit: 

^max ~ (3v/2ck) 1/2 qBAX ^ A(v/ck) 1/2 B a AX w TeV.(8) 

The corresponding dynamical ranges are 



^max/^min 



200(x/k) 1 / 2 (v/c) 1 ^ 2 B Q 3/2 AY 9 - 1 and (9) 



10(v/cfc) 1/2 AAio AYf 1 . 



(10) 



In general, to accelerate particles to TeV energies, at 
least mildly relativistic VA-values will be required. In addi- 
tion, for electrons, if converging flow acceleration is to be 
efficient, AYg < B^ 2 . This will typically require in as- 
trophysical sources either very high compression ratios, or 
very low magnetic fields, the latter implying a very large 
structure (e.g. magnetized turbulent ISM) if the process is 
to have observable radiative effects. For protons, the Hillas 
and diffusive limits lead to the more reasonable conditions 
AX > xAF and > (u/cfc) -1 ' 2 AY, respectively. Note nev- 
ertheless that electrons remaining within the reconnection 
layer may still be accelerated by electric fields generated in 
the reconnection process, resulting in energies < -Emm- In 
this case, the acce lerated energy distribution could also be 
Q(E) oc 1/E fe.g. lZenitani fc Hoshinoll200lD . 

Synchrotron emission from electrons injected with 
Q(E) oc 1/E and cooled through this process would have 
a spectral energy distribution v L v oc v 1 ! 2 . On the other 
hand, proton radiation would show v L v oc v for syn- 
chrotron (or oc v 1 / 2 in saturation) and for 7r°-decay from 
pp and photomeson production (above threshold). The last 
two processes also generate secondary pairs from n^- 
decay with spectra and luminosities similar to those of 
the photons from 7r°-decay. These pairs will predominantly 
cool through synchrotron, like primary electrons, also with 
v L v oc v 1 / 2 but peaking at much higher energies. 

Before exploring the case of magnetized jets, we note 
that our results may also be applied to weakly magne- 
tized flows with radiative shocks. Particle acceleration be- 
tween the upstream and the cooled downstream would re- 
quire sufficiently energetic particles to cross the adiabatic 
region of the post-shock flow and reach the cooled down- 
stream medium. Another condition would be that the parti- 
cle cooling/escape timescales were longer everywhere than 
the acceleration one. In this context, young stellar object 
jet shocks, or the dense winds of two massive stars colliding, 
may be sources worthy of being investigated. 

3. The case of magnetized jets 

Powerful astrophysical outflows, such as pulsar winds, 
and microquasar and AGN jets, are expected to be mag- 
netically d ominated at th e ir formation and acceleration 
zones fe.g. ICoronitil 119901: iBogovalov fc Tsinganosl 119991: 
iBeskin fc Nokhrinal 120061: iKomissarov et al.ll2007D . In pul- 
sar winds, magnetic reconnection may take place nat- 
urally in a curr ent sheet that originated at the win d 
equatorial region (|Coronitilll990l : ILvubarskv fc Kirkll200ll ). 
In relativistic jets, a -B-fie ld with reversals could come 
from the accretion disc (e.g. iBarkov fc Baushevll201lF) , but 
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magnetic dissipation may need to be triggered through 
jet acceleration, which could drive magnetic reconnec- 
tion through th e Kruskal-Schwarzschild instability (see 
Lvubarskvl 2010h . A similar though more extreme ef- 
fect can be expected if an obstacle is entrained by the 
iet (e .g. iHubbard k Blackmanl [200l lAraudo et all 1 2009 



20101: IBarkov et al.l l2010l: iPerucho k Bosch-Ramonl 12012 : 



Bosch-Ramon et al.|[20T2t iBarkov et alJl2012T T with a dom- 



inant B- field with polarity reversals of size AX. Under the 
ongoing magnetic reconnection, the incoming jet material 
and the obstacle itself can play the role of converging flows 
with a high AX/ AY ratio. To avoid the suppression of the 
acceleration process, the shocked obstacle must fulfil the 
following conditions: high inertia, to avoid quick dragging 
by the jet; not too high density, to avoid fast radiation cool- 
ing; and a mean free path much shorter than the obstacle 
size (-D ), to avoid fast particle escape. 

For a jet/obstacle interaction at Zj ~ 100i?s c h from the 
central object, where R Sch 3 x 10 13 (M/10 8 M Q ) cm, a jet 
width Dj ~ 0.1 Zj ~ 3 X 10 14 (M/10 8 M ) cm and magnetic 

field B ~ 1O 3 L J 1/ 4 2 4 (M/1O 8 M )- 1 G, and D Q > AX, the 
luminosity budget available for dissipation will be ~ 
(D Q / Dj) 2 Lj < Lj, where Lj is the jet power. Assuming 
AX ~ i?sch, the electron and proton maximum energies 
will be £&„ ~ l(M/lO 8 Af ) 1 / 2 L j r 4 1 4 /4 TeV, and ~ 

10 7 4«TbV. 



B 



Given the strong magnetic fields, the condition AYg < 
may not hold; if so, electrons could only be ac- 



-3/2 



1.5 



celerated within the reconnection layer. Thus, in AGN 
jet/obstacle interactions, magnetic reconnection itself could 
potentially accelerate electrons up to TeV energies. 
Protons, on the other hand, could be accelerated in the 
converging flow£| to ultra high e nergies or even beyond de- 
pending on AX, Lj, and D Q (see lGiannio 9 120101 for similar 
results). Electrons would yield hard synchrotron radiation 
peaking around 10 MeV. Proton synchrotron emission, of 
hard spectrum and peaking around 100 GeV , may be effi- 
cient as well under these conditions (see, e.g.. lBarkov et al.l 
l2012f ). Emi ssion through pp ins ide the obstacle could be also 
significant (|Barkov et al.ll2010h , as well as photomeson pro- 
duction in very bright AGN. Secondary very energetic 
pairs from 7r ± -decay would radiate ultra high-energy pho- 
tons through synchrotron emission that could pair-create 
in the jet radio fields on pc-scales, being reprocessed as 
synchrotron and IC photons of lower energies. 

In microquasar jets, where, say, M ~ 10 Mq and Lj ~ 
10 36 erg s _1 , magnetic reconnection in the jet/obstacle 
boundary at Zj ~ 100 i?sch could lead to synchrotron MeV 
flares produced by GeV electrons. A minor synchrotron self- 
Compton component in GeV may be also expected. PeV 
protons may be accelerated in the converging flows, and if 
they met a nearby thick target, for instance the obstacle it- 
self or a dense field of energetic photons, they could produce 
hard emission peaking at ~ 100 TeV through pp or pho- 
tomeson production, plus hard TeV synchrotron radiation 
from the subsequent pairs. In compact and powerful ob- 
jects photon-photon absorption would be severe, radiation 
being re-emitted as soft gamma rays via synchrotron emis- 
sion. It is not clear which kind of obstacle may be present in 
the innermost regions of a microquasar jet, although in the 



1 Even if jets were formed only by e pairs, protons may be 
entrained from the environment, e.g. the obstacle itself. 



(unlikely) event that the jet remained strongly magnetized 
up to the binary scales (3> 100 R^ch), a clumpy ste llar wind 
could provide those obstacles (|Araudo et al.l l2009). 

In both microquasars and AGN, the reconnection 
process can be efficient as long as the obstacle has a 
velocity different from that of the flow. In this pro- 
cess, for D a ~ Dj, a significant fraction of Lj may 
be released, with potential fluxes at high energies F ~ 
10~ 10 erg cm^s- 1 xL dM (d/lOOMpc)" 2 for AGN, and 
x ^d.36 (d/10 kpc) -2 for microquasars. This radiation would 
be roughly isotropic, but for obstacles small and light 
enough to be accelerated by the jet, the emission would be- 
come progressively Doppler-boosted and enhanced for spe- 
cific viewing angles. The obstacle may expand, enhancing 
the radiation, and possibly fragment while accelerated. All 
this should lead to complex sp ectral and variability patterns 
(e.g. IBarkov et al.ll2010l l2012f) . The typical timescale of the 
whole interaction proc ess will be the time the ob stacle re- 
mains as such (see, e.g. . iBosch- Ramon et al.l[2012h . whereas 
the reconnection events will have a variability timescale 
~ mm[D , AX]/v, where v < c. 
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